Quantum escape kinetics over a fluctuating barrier 
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The escape rate of a particle over a fluctuating barrier in a double well potential exhibits resonance 
at an optimum value of correlation time of fluctuation. This has been shown to be important in 
several variants of kinetic model of chemical reactions . We extend the analysis of this phenomenon 
of resonant activation to quantum domain to show how quantization significantly enhances resonant 
activation at low temperature due to tunneling. 



I. INTRODUCTION 



The enhancement of small periodic signal by noise in 
a nonhnear system has been a theme jl| of topical in- 
terest over more than two decades. Ever since the ob- 
servation of this phenomenon of stochastic resonance in 
varied theoretical and experimental contexts, a number 
of noise-induced resonance effects have been reported 
JJ- A prototypical effect of this kind, known 
as resonant activation discovered by Doering and Godua 
[l^ . concerns a resonance effect in the escape rate of a 
particle over a fluctuating barrier in a bistable potential, 
where the resonance [TElIlllllllllilliEllillll 
can be achieved by observing the variation of mean first 
passage time as function of flipping rate of fluctuation 
of the barrier height. In the simplest possible term the 
phenomenon can be realized in a model with a linear 
barrier with a slope which fluctuates between two values. 
The phenomenon has triggered a lot of theoretical activ- 
ity around the Markovian and non-Markovian variants of 
kinetic models for chemical reactions 0, 0, 0, 0| and 
several related issues and been experimentally observed 
by Mantegna and Spagnolo 22]. This fluctuation of po- 
tential is also important in various problems of chemical 
physics, notably in molecular dissociation dynamics |23j |. 
protein folding 24J among others. The purpose of this 
paper is to extend the analysis of resonant activation to 
quantum domain. Since at low temperature thermal ac- 
tivation is accompanied by tunneling, the question that 
naturally arises is how noise-induced resonance effects 
manifest themselves in a quantum system. Specifically 
our object here is two-fold: First, we intend to under- 
stand the counterintuitive role of external noise in bring- 
ing out the resonance behaviour in a bistable quantum 
system where the mean escape time is varied as a func- 
tion of correlation time of noise. Second, it is worthwhile 
to examine the nature of this resonance in presence of 
generic quantum effects like tunneling at low temperature 
and allow ourselves a fair comparison with the classical 
results. In what follows we carry out a theoretical and 
a numerical study of quantum stochastic dynamics of a 
bistable system with a barrier height fluctuating due to 
an external Ornstein-Uhlenbeck noise and show how the 



mean escape time is profoundly influenced by the statis- 
tical properties of the noise particularly the correlation 
time in exhibiting the quantum resonance activation. 



II. QUANTUM STOCHASTIC DYNAMICS 

A. General aspects 

To derive quantum Langevin equation from a mi- 
croscopic picture consider the well-known standard 
system-reservoir model with following form of the 
HamiltonianlSSl 



N 

E 



(2.1) 



Here x and p are the coordinate and momentum op- 
erators of the particle and {qj,Pj} are the set of coor- 
dinate and momentum operators for the reservoir oscil- 
lators coupled linearly through the coupling constants 
Hj[j — 1, 2, ...). The potential V{x, t) is due to the exter- 
nal force field for the Brownian particle. The coordinate 
and momentum operators follow the usual commutation 
rules {x,p\ = ih and {qi,pj} = ihSij. Eliminating the 
bath degrees of freedom in the usual way we obtain the 
operator Langevin equation for the particle 



mS:+ dt'-f{t~t')S:{t') + V'{x,t) ^ F{t) (2.2) 



where the noise operator F(t) and the memory kernel 
7(t) are given by 



F(t) = ^ [^{9j(0) — x{0)}Kj cosLUjt + Ky^pj{0) sinujji 



and 



f{t) — Kj COS UJjt 



(2.3) 



(2.4) 
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respectively with kj — uj'j. On the basis of quantum 
mechanical average (...) over the bath modes with coher- 
ent states and over the system mode with an arbitrary 
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state Eq.(2.2) can be cast into the form of the generahzed so that the statistical averages {■■■)s over the quantum 
quantum Laneevin equation [H IH El US IM HI HI mechanical mean value Oj of the bath variables are de- 
m- fined as 



mx+ f dt'-f{t-t')x{t') + V'{x,t) = f{t) + Q{x,{Sx'')) 
Jo 

(2.5) 

where the quantum mechanical mean value of the po- 
sition operator (x) = x. Here the quantum dispersion 
term Q to the potential, is given by 



Qix, {Sxn)^V'ix,t)-{V'ix,t)) 



(2.6) 



which by expressing x{t) — x{t) + 5x{t) in V{x, t) and 
using a Taylor series expansion around x may be rewrit- 
ten as 



Q{x,{5x-)) = -Y,-^V^^+'\x,t){5x^ 



(2.7) 



n>2 



Here the (n -I- l)th derivative with respect to 

X. The calculation of Q rests on the quantum correc- 
tion terms ((Jx") which can be calculated order by order 
by solving a set of quantum correction cquations(as dis- 
cussed in the later part of this section) . Furthermore the 
quantum mechanical mean Langevin force is given by 



(O,). = j O, P, d{p,{Q)) d{{q,m - (x(0))} (2.12) 

Here n{uj) is given by Bose-Einstein distributions 
(glr — i)-i. is the exact solution of Wigner equa- 
tion for harmonic oscillator [stL [s^ and forms the basis 
for description of the quantum noise characteristics of the 
bath kept in thermal equilibrium at temperature T. N is 
the normalization constant. In the continuum limit the 
fluctuation-dissipation relation (|2.1U|) can be written as 



(/W/(i')> 
1 



huj 



dco k(uj) p{uj) Two coth(— — ) cosa;(i — t') 



2kT' 



(2.13) 



where we have introduced the density of the modes 
p{uj). Since we are interested in the Markovian limit 
in the present context, we assume K,{u})p{u)) — -7, and 
Eq.(EiEJl then yields 



{f{t)f{t'))s = 2D,5{t^t') 



(2.14) 



f{t) — ^ |^((7j(0)) — (x(0))Kj cosLUjt + Ky'^pj{0) smu jt 

(2.8) 

which must satisfy noise characteristics of the bath at 
equilibrium , 



with 



{f{t))s 
{f(.t)f{t'))s 



(2.9) 
'^"^ huji I coth cosijj, (t — t') 

(2.10) 



j 



Ea. H2.10() expresses the quantum fluctuation- 
dissipation relation. The above conditions Ea. H2.9|l - 
Ea. (|2.1U|l can be fulfilled provided the initial shifted 
co-ordinates {{qj{0)) — (x(0))} and momenta {pj{0)) 
of the bath oscillators are distributed according to the 
canonical thermal Wigner distribution jS^, l38| of the 
form 



Dn — —"/hun coth — 

^ 2' " 2kT 



(2.15) 



(The passage from Ea. (|2.13|) to Ea. (|2.14|l is given in the 
appendix A) 

0^0 refers to the static frequency limit. Furthermore 
from Ea. (|2.4() in the continuum limit we have 



^{t-t') = ^5{t-t') 



(2.16) 



7 is the dissipation constant in the Markovian limit. 
In this hmit Ea. (|2.5|) therefore reduces to 



mi + -ix + V'{x, t) = f{t) + Q{x, ((5£")) (2.17) 

In order to consider the activated processes in a bistable 
potential with a fluctuating barrier height we consider 
the potential of the form 



PMm)-{m)].{pM)) 
A^exp^ WM? + i^3mm~{m)]' 



(2.11) 



V{x,t) = U{x)+g{x)at), 



(2.18) 



where U{x) is a bistable potential (- 



a^2 , 6^4 



x , a, b being 



constants ) with a barrier at the metastable point a; = 
and two stable points at x = ±(|-)2 . The fluctuations in 
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the potential are driven by an Ornstein-Uhlenbeck noise 
process 



(2.19) 



where r]{t) is zero mean (5-correlated Gaussian noise. The 
stochastic process £,{t), can be characterized by the fol- 
lowing set of equations. The probability function /o(^), 
the variance and the correlation function of the noise 
are given 



p(0 = yWexp(-^) 



where 



and 



/+00 
-oo 



(2.20) 



(2.21) 



(2.22) 



respectively. The barrier can also be subjected to di- 
chotomic fluctuation ^(t) = {—a,P} that flips be- 
tween two values with flipping rate fia and /i^ respec- 
tively. This process can also be taylored as zero mean 
valued, exponentially correlated and with a Gaussian dis- 
tribution. 



B. Overdamped limit 

To proceed further we now confine ourselves to over- 
damped condition so that the quantum Langevin equa- 
tion in one variable takes the following form: 

-fx^ax- bx^ + g'ix)^{t) + fit) + Q{x, ((5£")) (2.23) 

The above classical-like stochastic differential equation 
contains the essential quantum features, through the 
terms flf) which represent the quantum noise of the heat 
bath and the another term Q{x, {Sx^)) which essentially 
arises due to nonlinear part of the system potential. The 
nonlinearity and the quantum effects are entangled in the 
latter quantity modifying the classical part of the poten- 
tial. Thus the classical potential force —V'{x^t) is mod- 
ified by the quantum dispersion. In absence of quantum 
dispersion term Q{x, ((5a;")) and with Dq — > ^kT as one 
approaches the classical limit {kT 3> ^o); the quantum 
Langevin equation reduces to classical one. A proba- 
bilistic description of the system Ea. (|2.23|l is provided 
by the time-dependent probability distribution function 
P{x, ^, t). For a finite correlation time of the noise driving 
the barrier to fast fluctuations it is possible |23| to go over 
to an approximate (exact in the limit r — > 0) description 
which allows one to transform externally driven Langevin 



equation Ea. H2.23() into an autonomous Langevin equa- 
tion for the position variable x. The corresponding equa- 
tion for probability density function P{x, t) is given by 



dP{x,t) 
di 



dx 

+ a^ — 

dx 



[n{x,T)f{x)]P{x,t) 
d 

k{x, T)g'{x) — K{x, T)g'{x) 



k{x,t) — k{x,t) 



Pix,t) 
P{x,t) (2.24) 



where 



1 



assumed 
1 

and 



to unity, k{x, t) = 
fix) = -U'ix) + 



Qix,{Sx"')). The quantum nature of the above 
Fokker-Planck equation governing evolution of Pix,t) 
is manifested through the quantum correction to the 
potential term and the quantum diffusion coefficient 
characterizing the thermal bath. 

The quantity of special interest here is the mean ther- 
mally activated escape time of the particle from the well 
whose barrier is subjected to fluctuation. The particle is 
governed by the stochastic dynamics Eq. H2.24() . It is im- 
portant to emphasize that the typical mean escape time 
should be much larger than the time scale of the deter- 
ministic quantum dynamics i — ~U'ix) + Qix, {Sx'^)) for 
nonzero fluctuation of the potential well, and it is further 
required that the strengths of thermal noise and barrier 
fluctuation must be small in comparison with the barrier 
height. The solution of Ea. (|2.24|l in the stationary state 
reads as 



Psix,T) = 
I 



■g'ix)'- 



X exp 



/(y) 



1 - 



V 



DM^^g'ivY 



dy 



(2.25) 



The stationary probability distribution function Ps(x,t) 
of the stochastic process x(t), which is the quantum 
mechanical mean position, is bimodal in nature with 
nonzero probability current at the barrier top even at 
zero temperature characterizing the zero point contri- 
bution of the thermal bath (as T ^ 0, Dg hujo). 
Pix, t) is additionally modified over its classical nature 
by a quantum contribution characterizing a correction 
due to anharmonic part of the system potential. 

To proceed further it is necessary to find out the quan- 
tum correction term Qix, {dx")) more explicitly. To this 
end we return to the overdamped operator equation H2.2|) 
and use xit) = xit) + Sxit), where a;(i)(= (i(t))) is the 
quantum mechanical mean value of the operator x. By 
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construction [5x, 5p] = ih and (Sx) = 0. We then ob- 
tain the quantum correction equation in the overdamped 
hmit after quantum mechanical averaging with the co- 
herent states over the bath operators as 



V"(x,t)Sx + y ^y"+i(a;,i)((5x" - (Sx'')) = 



With the help of operator equation 
the equations for ((5x") 



(2.26) 
we obtain 



1 



dt 



(Sx^) = — [2V"{x, t){Sx^) + V"'ix, t){Sx'^)] (2.27) 



W"ix,t){Sx'^) + -V"'{x,t){Sx^) 



±.V"\x,t){6xr 



(2.28) 



and so on. To take into account of the leading order 
contribution {6x'^) explicitly we may write 

d{Sx^) = ~-V"ix, t){5x^)dt (2.29) 
7 

The overdamped deterministic motion on the other hand 
gives ^dx = —V'{x,t)dt which when used in Ea. H2.29|) 
yields after integration 

{6x^) ^ A,[V'ix,t)f (2.30) 
Aq is the quantum correction parameter, as given by 
Aq — |^77^"^7)p- 1 being a given quantum mechanical 
mean position. 



III. RESONANT ACTIVATION IN THE 
QUANTUM SYSTEM 

We are now in a position to analyze the quantum res- 
onant activation. In our present problem the fluctuating 
part of the potential associated with harmonic term is 
assumed to be of the form g{x) = ^ and x 



±,/f and 

b 



X = are the absolute minima and maximum of the po- 
tential U(x) respectively. Having known the stationary 
distribution Eq. (|2.25|) along with the quantum correction 
in f{x) the calculation of stochastic dynamics is straight- 
forward. This can be obtained using the standard result 
on mean first passage time .26i] as 



(T) 



dx 



D{x,t)Ps{x) J_ 



where Dejfix, t) = Dq {l + Rg'{x)^) 
to give 



Ps{x)dy (3.1) 



27^1^^^ 
aV2(l-A2) 



AT4//(i?,r) 



D„ 



(3.2) 



Here Ai and A2 both are the leading order quantum 
correction terms as given by 



Ai 



f Q{x,{Sx^)) 
V 5'(^) 



(3.3) 



(3.4) 



and R—jy\ AVeff{R,T) is the effective barrier height 
of the following form 

AK//(i?,r) 



{-U'{x) + Q) 


i+v(.)(-^); 


dx 




l + Rg'{x)^) 



(3.5) 

Ea. (|3.2() is the central analytical result of this paper. To 
analyse the theoretical results on the essential features 
of the activated escape following quantum stochastic dy- 
namics over a fluctuating potential well we now resort 
to numerical simulation of the Eas. (|2.23|l . (|2.19() and the 
quantum correction equations H2.27|l and H2.28II simulta- 
neously using standard Heun's algorithom. A very small 
time step (At) 0.001 for numerical integration has been 
used. In our simulation we follow the dynamics of each 
particle starting in the left well at a; = — till it arrives 
at the barrier top at a; = 0. The first passage time being 
a statistical quantity due to the random force we calcu- 
late the statistical average of the first passage time over 
1,000 trajectories. We present the numerical results in 
Fig.l to Fig. 4 for different parameters such as tempera- 
tures, barrier height of the potential and noise strength of 
the Ornstein-Unhenbeck noise process. All the curves ex- 
hibit minimum at optimal r values. Increasing strength 
of thermal or non-thermal noise results in enhancement 
of escape rate and the minima are shifted towards the ori- 
gin. Physically this implies that with increase in Dg or 
cr^, the escape time over the barrier changes as a result 
of which the switching time of the barrier fluctuations 
matches the order of escape time at lower r values. In 
order to examine the influence of the barrier height on 
resonant activation we plot in Fig. 3 the variation of mean 
escape time as function of r for several values of barrier 
heights. Increase in the barrier height results in enhance- 
ment of escape time and the minima are shifted to higher 
T values. 

We find that the barrier height Eq. H3.5|l monotonically 
decreases to a limiting value with increasing t. Due to 
~^the presence of system nonlinearity the barrier height is 
slightly modified by an added slope. It is thus expected 
that during the temporary stay of the particle close to a 
minimum it makes unsuccessful attempts to escape and 
once the escape takes place it should occur at a slower 
rate in quantum case due to added slope to the effective 
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potential than the corresponding classical case. However, 
our theoretical and numerical calculations reveal that the 
mean escape time is lower in magnitude for the quantum 
particle at low temperature and the difference becomes 
insignificant at higher temperature. This is shown in the 
Fig. 4. This behaviour can be interpreted in terms of an 
interplay between the quantum diffusion coefficient Dq 
and the quantum correction due to system nonlinearity 
appearing in Vef f{R,T). When the temperature of the 
system is very low, i. e., in the vacuum limit or in the 
deep tunneling region the anharmonic terms in the po- 
tential do not contribute significantly. On the other hand 
as temperature of the system increases significantly, Dq 
increases resulting in decrease of the effective potential 
and hence Dq and Q compete to cancel the effect of each 
other at higher temperature. Finally in Fig. 5 we make 
a comparison between the numerical simulation and the 
corresponding theoretical result. The theoretical result 
agrees fairly well with our numerical result. 

In addition to the barrier height the prefactor is also 
affected by the quantum correction term. The approxi- 
mate theoretical mean escape time Ea. (|3.2|l implies that 
both the barrier height and the frequency factor have an 
important role to play with the activated escape process. 
We are now in a position to point out that the ther- 
mally activated resonance phenomenon is controlled by 
three different time scales as emphasized earlier . The 
three relavent time scales are the corresponding time of 
barrier fluctuation (0 <C r <C oo), typical escape time T, 
which is much larger than the time scale describing de- 
terministic motion i = —U'{x) -I- Q'{x, ((5x")), the third 
time scale being Ta of the escape attempt {Ta ^ T). 



IV. CONCLUSION 

Based on the theoretical study and numerical simu- 
lation of quantum stochastic dynamics in a double-well 
system with a fluctuating barrier under influence of a 
Gaussian color noise, we have examined the phenomenon 
of resonant activation. The governing equations are clas- 
sical looking in form but quantum mechanical in their 
content. A key point of the present analysis is to describe 
the thermal bath in terms of a canonical thermal Wigner 
distribution of harmonic oscillators. This distribution re- 
main positive definite even at absolute zero signifying a 
pure state and allows us to look for the external noise- 
induced resonance when the generic quantum effects in 
the system make their presence felt. We have shown that 
quantization significantly enhances resonant activation at 
low temperature due to tunneling. For higher tempera- 
ture as well as for stronger noise strength the resonance 
effects get more pronounced. Since tunneling accompa- 
nies activation, it is expected that the resonant activation 
can be observed even at absolute zero, where the reso- 
nance effect essentially is due to the vacuum field. 



APPENDIX A: (THE PAS SAGE FROM EO. (E351l 
TO EQ. i!S.14ll ') 

We start from basic definition(Louisell p-426) js^ 



i+At pt+At 

dt 



dt' (/W/(t')>. (Al) 



Using Eg. lfTT^ in Eq.||JTJ yields 



- 2At 



— [ duj k(uj) p{uj)hLj coth ( 

V2fcr 

t+At pt+At 

XI dt I dt' cosuj{t - t') 



(A2) 



Explicit integration over time gives 



2D„ 



where 



f 

2Ai 



' dw k{ll}) p{uj)fiw coth — — At) 
\2kT J 

(A3) 



4 2 

sm — — 



(A4) 



Putting k{uj) p{lo) = 7, we obtain 



2Dq = 



(- 



huj \ sin 



2 u:At 



do. ^u. coth j (A5) 



Following Louisell(p —426) |33 we have under Marko- 
vian condition, the correlation time Tc ^ At, the 
coarse-grain time(over which the probability distribution 
function evolves) 

Thus as At oo(in scale of Tc which goes to zero) the 

sin^ 

function oscillates violently so that one takes the 

V 2 > 

slowly varying quantity [hw coth out of the integra- 
tion over frequency with an average value fuvo coth , 
ljq be an average static frequency. Since the inte- 
gral """^a^^ rfx = TrAt it follows immediately from 
Ea.(lA5^ 



2Dq — 7 hhjQ coth 



V2fcr 



(A6) 



as given in Ea. (|2.15|) . 

Again starting from Ea. H2.13|l . we use the same argu- 
ment as before to have 



{f{t)f(t'))s 

f { huj \ 

= — / du! K{uj)p{uj)huj coth J cosaj(t — t') 

2 Jq \2k1 J 

and we use 
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duj cos UJT — TT S{t) 



{f{t)f{t'))s = 2DgS{t-t') 



to obtain 



{fit).fit')). 
1 





'2 




1 doj 


-7 

TT 


nu! cotn — — 


h 


\2kT J 



2 

7 f^o 
coth 

TT 

7 hujo coth 



(SI)'*"-''' 



(A7) 



Therefore from Ea. (|A6|l and Ea. (|A7|) we have 



which is Ea. (|2.14|l . 

Thus the derivation within Markovian approximation 
clearly depends on the time scale separation. The results 
are valid even at absolute zero as emphasized by Louisell. 



Acknowledgments 

Thanks are due to the Council of Scientific and in- 
dustrial research, Govt, of India, for partial financial 
support. 



[1] R. Benzi, G. Parisi, A. Sutera, and A. Vulpiani, Tellus 

34, 10 (1982); R. Benzi, A. Sutera, and A. Vulpiani, J. 

Phys. A 14 L 453 (1981) 
[2] B. McNamara and K. Wiesenfeld, Phys. Rev. A, 39, 4854 

(1988); L. Gammaitoni, F. Marchesoni, E. Menichella- 

Saetta and S. Santucci, Phys. Rev. Lett. 62, 349 (1989). 
[3] B. McNamara, K. Wiesenfeld and R. Roy, Phys. Rev. 

Lett. 60, 2626 (1988) 
[4] L. Gammaitoni, P. Hanggi, P. Jung and F. Marchesoni, 

Rev. Mod. Phys., 70, 223(1998) 
[5] F. Julicher, A. Ajdari and J. Prost, Rev. Mod. Phys. 69, 

1269 (1997). 

[6] M. O. Magnasco, Phys. Rev. Lett. 71, 1477 (1993), C. 

R. Doering, W. Horsthemke and J. Riordan, Phys. Rev. 

Lett. 72, 2984 (1994). 
[7] H. Linke, T. E. Humphrey, A. Lfgren, A. O. Sushkov, 

R. Newbury, R. P. Taylor, and P. Omling, Science 286, 

2314 (1999); L. Q. Zhou, X. Jia and Q. Ouyang, Phys. 

Rev. Lett. 88, 138301(2002) 
[8] J. Garcia-Ojalvo, A. Hernandez-Machado and J. M. San- 

cho, Phys. Rev. Lett. 71, 1542 (1993) 
[9] A. Becker and L. Kramer, Phys. Rev. Lett. 73, 955 

(1994), A. Becker and L. Kramer, Physica D 90, 

408(1995) 

[10] J. M. R. Parrondo, C. Van den Broeck, J. Buceta and F. 

J. de la Rubia, Physica A 224, 153(1996) 
[11] A. A. Zaikin and L. Schimansky-Geier, Phys. Rev. E. 58, 

4355(1998) 

[12] J. Buceta, M. Ibanes, J. M. Sancho and K. Lindenberg, 

Phys. Rev. E 67, 021113(2003) 
[13] S. Dutta, S. S. Riaz and D. S. Ray, Phys. Rev. E 71, 

036216 (2005) 

[14] C. R. Doering and J. C. Gadoua, Phys. Rev. Lett. 69, 
2318 (1992) 

[15] M. Bier and R. D. Astumian, Phys. Rev. Lett. 71, 1649 

(1993) ; U. Ziircher and C. R. Doering, Phys. Rev. E 47, 
3862 (1993) 

[16] C. Van den Broeck, Phys. Rev. E 47, 4579 (1994) 

[17] J. J. Brey and J. Casado-Pascual, Phys. Rev. E. 50, 116 

(1994) ; O. Flomenbom and J. Klafter, Phys. Rev. E. 69, 



051109 (2004) 

[18] T. D. Shepherd and R. Hernandez J. Chem. Phys. 115, 

2430 (2001); A. M. Berezhkovskii, A. Szabo, G. H. Weiss, 

and H. Zhou J. Chem. Phys. Ill, 9952 (1999) 
[19] M. Boquria, J. M. Porra, J. Masoliver and K. Lindenberg, 

Phys. Rev. E 57, 003990 (1998) 
[20] M. Marchi, F. Marchesoni, L. Gammaitoni, E. Menichella 

-Saetta, and S. Santucci Phys. Rev. E 54, 3479 (1995) 
[21] P. Hanggi, Chem. Phys., 180 157 (1994); R. F. Fox, Phys. 

Rev. A. 34, R4525 (1986); Y. Jia and J. R. Li Phys. Rev. 

E. 53, 5786 (1996); X. Luo and S. Zhu, Phys. Rev. E 67, 

021104 (2003) 

[22] R. N. Mantegna and B. Spagnolo, Phys. Rev. Lett. 84, 
3025 (2000) 

[23] J. Maddox , Nature(London) 359, 771 (1992). 

[24] J. Wang and P. Wolynes, Comput. Phys. 180 141 (1994) 

[25] J. M. Martinis, M. H. Devoret and J. Clarke, Phys. Rev. 

Lett. 55, 1543 (1985) 
[26] R. L. Straonovich, Topics in the Theory of Random 

Noise, Vol. 1 (Gordon and Breach, New York, 1963); K. 

Lindenberg and B.J. West, J. Stat. Phys. 42, 201 (1986); 

J. Masoliver, B.J. West, and K. Lindenberg, Phys. Rev. 

A 35, 3086 (1987); D. Mei, G. Xie, Li Cao, and D. Wu, 

Phys. Rev. E 59, 3880 (1999) 
[27] D. Banerjee, B. C. Bag, S. K. Banik and D. S. Ray, Phys. 

Rev. E 65, 021109 (2002). 
[28] D. Barik and D. S. Ray, J. Chem. Phys. 119, 12973 

(2003). 

[29] D. Barik, S. K. Banik and D. S. Ray, J. Chem. Phys. 

119, 680 (2003). 
[30] D. Banerjee, B. C. Bag, S. K. Banik and D. S. Ray, J. 

Chem. Phys. 120, 8960 (2004). 
[31] D. Barik and D. S. Ray J. Chem. Phys. 121, 1681 (2004). 
[32] P. K. Ghosh, D. Barik and D. S. Ray, Phys. Rev. E 71, 

041107 (2005) 

[33] P. K. Ghosh, D. Barik and D. S. Ray, Physics Letters A 
342 12 (2005). 

[34] D. Barik and D. S. Ray, J. Stat. Phys. (2005) To be 
published. 

[35] R. Zwanzig, J. Stat. Phys.9, 215 (1973). 



7 



[36] A. Mielke, Ann. Phys. 4, 476 (1995); J. M. Sancho, Phys. 
Rev. A 31, R3523 (1985); J. M. Gutirrez, A. Iglesias, and 
M. A. Rodrgucz, Phys. Rev. E. 48, 2507 (1993) 

[37] E. P. Wigner, Phys. Rev. 40, 749 (1932). 

[38] M. HiUery, R. F. O'Connell, M. O. Scully and E. P. 
Wigner, Phys. Reps. 106, 121 (1984). 

[39] W. H. Louisell, Quantum Statistical Properties of Radia- 
tion, (John Wiley and Sons, New York, 1990); 

Figure Captions 

Fig.l. A plot of mean escape time vs t using Gaussian 
colour noise for different temperatures [(i) T = 0.1 (dot- 
ted line), (ii) T = 1.0 (dashed line), (iii) T = 2.0 (solid 
line)] and the other parameter set a = 0.5, b = 0.005 and 
a2 = 1.0. 

Fig. 2. A plot of mean escape time vs t using Gaussian 
colour noise for different cr^ [(i)f^ = 0.5 (dashed line), 

(ii) CT^ = 1.0 (dotted line), (iii) ct^ = 1.5 (solid line)] and 
the other parameter set a = 0.5, b = 0.005 and T = 1.0. 



Fig. 3. A plot of mean escape time vs r using Gaus- 
sian colour noise for different barrier heights [(i) a = 0.3 
(dashed line), (ii) a = 0.5 (dotted line), (iii) a = 0.7 (solid 
line)] and the other parameter set b = 0.005, = 1.0 
and T = 1.0. 

Fig. 4. A comparison between classical(dotted line) and 
quantum (solid line) mean escape time describing reso- 
nant activation phenomenon using Gaussian colour noise 
for the parameter set a = 0.5, b = 0.005 and cr^ = 1.0 at 
two different temperatures T = 1.5 and T = 0.1. 

Fig. 5. A comparison between numerical (dashed dot 
dot line and dashed line) and analytical (solid and dotted 
line) results using Gaussian colour noise for different 
[(i) = 0.5) (solid and dashed dot dot line), (ii) = 1.5 
(dotted and dashed line)] and the other parameter set 
a = 0.5, b = 0.005 and T = 1.5. 
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